arXiv: 1509.00119vl [physics.optics] 1 Sep 2015 


Elusive pure anapole excitation in homogenous spherical 
nanoparticles with radial anisotropy 

Wei Liu/’* Bing Lei/ Jianhua Shi/ Haojun Hu/ and Andrey E. Miroshnichenko^ 

^ College of Optoelectronic Science and Engineering, National University of Defense Technology, Changsha, 

Hunan 410073, P. R. China. 

^Nonlinear Physics Centre, Research School of Physics and Engineering, Australian National University, 

Acton, ACT 0200, Australia. 

* Corresponding author: wei.liu.pku@gmail.com 


Copyright © 2015 Wei Liu et al. This is an open access article distributed under the 
Creative Commons Attribution License, which permits unrestricted use, distribution, and 
reproduction in any medium, provided the original work is properly cited. 


For homogenous isotropic dielectric nanospheres with incident plane waves, Cartesian 
electric and toroidal dipoles can be tunned to cancel each other in terms of far-held 
scattering, leading to the effective anopole excitation. At the same time however, other 
multipoles such as magnetic dipoles with comparable scattered power are simultanesouly 
excited, mixing with the anopole and leading to a non-negligible total scattering cross 
section. Here we show that for homogenous dielectric nanospheres, radial anisotropy can be 
employed to signihcantly suppress the other multipole excitation, which at the same time 
does not compromise the property of complete scattering cancallation between Cartesian 
electric and toroidal dipoles. This enables an elusive pure anopole excitation within radially 
anisotropic dielectric nanospheres, which may shed new light to many scattering related 
fundamental researches and applications. 

1. Introduction 

Since the discovery by J. J. Thomson that electrons are constituents of atoms [1], there had been a lot 
of searches for charge-current distributions in which no electromagnetic energy are radiated, which can be 
applied to explain the stability of atoms [2]. This is due to the existing dogma since the establishment 
of Maxwell equations that accelerating charges would lose energy continuously through electromagnetic 
radiation, and consequently to prevent the collapse of atoms, the charges within them must be arranged 
into specific configurations. Such special charge-current configurations are also termed as non-radiating 
sources [2,3]. Those non-radiating sources are deeply related to the concepts of anapoles in quantum field 
theory, which have their roots in the law of parity unconservation [4], and can be used as possible explanations 
for the existence and features of neutrinos and cold dark matter [5,6]. 

For time-varying oscillating charge-current distributions, the dynamic version of anapoles is highly related 
to the so called toroidal multipoles, which recently have attracted enormous attention [7-20]. As indepen¬ 
dent terms in the multipole expansion of electromagnetic fields besides conventional cartesian magnetic 
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and electric multipoles [3], toroidal multipoles correspond to significantly different near-field charge-current 
distributions and will inevitably become significant when the operating spectral regime goes beyond the elec¬ 
trostatic region [20-22]. In the far field however, the toroidal multipoles have identical scattering patterns 
compared with their conventional Cartesian electric and magnetic counterparts, and that is to say, they are 
indistinguishable in terms of far-field scattering [7,13,21,23]. This at the same time also means that when 
properly engineered, the contributions of conventional and toroidal multipoles to the scattered field can be 
tuned to the same amplitude and out of phase, which leads to the far-field scattering cancellation and thus 
generate the so called electrodynamic anapoles [11,15-17,19,21,24]. The simplest case of such anapoles 
is induced by the destructive cancellation of Cartesian electric dipole and toroidal dipole, which has been 
demonstrated firstly in composite carefully engineered structures consisting of particle clusters [11,15]. Re¬ 
cently it has also been theoretically proposed and experimentally demonstrated in surprise that the anapole 
can be effectively excited even in individual homogenous dielectric nanoparticles [16]. However, the problem 
of the above demonstration in simple structures is that besides the anapole excitation, other multipoles 
(such as magnetic dipoles) of similar magnitudes in terms of scattered power are also excited and mixed with 
the excited electric and toroidal dipoles [16,17]. Such problems can be conquered within single core-shell 
nanoparticles of fundamental shapes (spherical or cylindrical) where the extra multipole excitations can be 
significantly suppressed, enabling elusive pure anapole excitations [17,19]. 

Here we propose an alternative approach and demonstrate that elusive pure anapole can be observed in 
homogenous radially anisotropic dielectric spherical nanoparticles under incident plane waves. It is shown 
that when the radial anisotropy of material permittivity is employed, the Cartesian electric and toroidal 
dipoles can be tuned to the spectral positions where other multipoles are significantly suppressed, without 
compromising the feature of complete destructive interference between them. Since Cartesian electric and 
toroidal dipoles can totally cancel each other in terms of far-field scattering and at the same time the 
excitation of other multipoles is negligible, it is proved that effectively pure anapole has been obtained. 
The anapole excitation is highly related to the exotic scattering phenomena of invisibility and possible 
applications of non-invasive sensing and detections, which may pave the way for many further scattering 
related fundamental researches and applications. 

2. Methods and Expressions 

A. Multipole Expansions with toroidal components 

For any charge-current distribution, when the geometric size of the source is far smaller than the effective 
wavelength (the retardation effect within the source can be neglected), the quasi-static approximation can 
be applied and then the radiated fields of the source can be expanded into two sets of conventional Cartesian 
electric and magnetic multipoles [3,21]. The lowest order terms of the two sets of multipoles are electric and 
magnetic dipoles, which can be expressed respectively as [throughout the paper we adopt the exp(—iwt-l-ik-r) 
notation for electromagnetic waves]: 

P=- f d^r3(r), M=— [ d^r[r xJ(r)], (1) 

-tu} J 2c J 

where c is the speed of light, uj is the angular frequency and J is the current. When the retardation effect 
within the source cannot be neglected (the size of the source gets larger and/or the effective wavelength 
gets smaller), the two sets of multipoles are not complete and thus other multipoles have to be included as 
correcting terms [20-22]. The lowest order correcting multipoles are toroidal multipoles, with the toroidal 
dipole expressed as as [7,21,23,24]: 

T = ^J d^r[{r-3{r))r-2r^J]. (2) 
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The radiated power of P, T and M are respectively: 
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(3) 

(4) 

(5) 


where /xq is the vacuum permeability and k is the angular wave number. We emphasize that the expressions 
above are for individual multipoles only and a simple sum of them does not lead to the total scattered power 
as there is interference between them [7,21,23] [see also Eq. (10) below]. It is also worth mentioning that 
the incorporation of toroidal multipoles does not make the expansion complete [20-22]. When the operating 
spectral region is far beyond the quasi-static regime, besides toroidal multipoles, other higher order correcting 
multipoles have to be also included [20-22]. 

For scattering particles, based on the charge conservation relation, the current (displacement current) 
distribution is 

J(r) =-iweo[e(r) - l]E(r), (6) 


where eo is the vacuum permittivity; £(r) is the relative permittivity, which is a scalar for isotropic particles 
or tensor for anisotropic particles. With the current distribution shown in Eq.(6), all the Cartesian multipoles 
including P, T and M and their scattered power can be directly calculated [see Eqs.(l)-(5)]. 


B. Electric and magnetic multipoles deduced from vector spherical harmonics of Mie scattering particles 

For spherical (both isotropic and radially anisotropic) particles with incident plane waves of electric field Eq 
and magnetic field Hq, both the fields inside the particle and the scattered fields can be expressed in transverse 
electric and transverse magnetic (in terms of electric field distribution)vector spherical harmonics [20,25,25- 
28]. On one hand, according to Eqs.(l)-(6), all the Cartesian multipoles including P, T and M and their 
scattered power can be analytically expressed and numerically calculated. On the other hand, based on the 
far-held equivalence, the scattered Helds can be viewed as the radiated Helds of a combination of conventional 
Cartesian electric and magnetic multipoles (which correspond to transverse magnetic and transverse electric 
vector spherical harmonics respectively). As a result the far-Held deduced electric dipole P(ai) and magnetic 
dipole M(6i) can be expressed respectively as [16,17,20,25,28-32]: 


, 67riai^ , diribi^^ 

P(ai) = Eq——Eq, M(6i) =/xq— 


(7) 


where ai and bi are the Mie scattering coefficients of the lowest order. The total scattered power is [25-28]: 


I ^ 12 oo 

= V(2m+l)(|a^|V|5^|'). (8) 

where Om and b^ are the Mie scattering coefficients of order m and Em = i'^E^ m{^+i) ■ clear that the 
scattered power of P(ai) and M(6i) is respectively: 
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C. Anapole excitation induced by destructive interference between P and T 

As is discussed above, P(ai) shown in Eq.( 7) has been obtained based on the complete vector spherical 
harmonic expansion and thus has included all the transverse magnetic dipolar scattering components [16,17, 
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Fig. 1: Schematic illustration of anopole excitation in a homogenous particle with plane wave incidence. The radius 
of the particle is R, and the particle is isotropic with refractive index n, or is radially anisotropic with transverse 
permittivity £t = and radial permittivity £r (the anisotropy parameter is expressed as: t] = etjer)- Within the 
particle, both T and P are excited, out of phase and of the same magnitude with regard to the scattered power. At 
the same time, the typical current J and magnetic held H distributions for T only have been shown. 


20,28-32]. While the conventional Cartesian electric dipole P shown in Eq. (1) is calculated through current 
integration and is only the lowest order expansion of P(ai) [20-22]. For operating spectral region beyond the 
quasi-static regime, the next order correcting term, which is the toroidal dipole T shown in Eq.( 7), should 
also be included. When higher order correcting terms beyond toroidal terms can be neglected, the scattered 
dipolar transverse magnetic field can be expressed as as [16,17,21,23]: 


E 


rad — 


exp(i/cr) 
47rreo 


n X P(ai) X n 


exp(ffcr) 

Airreo 


n X (P -b ikT) X n. 


( 10 ) 


where n is the unit vector along r. The existence of the anapole lies in P(ai) « P -I- ikT = 0 which can be 
satisfied when: 

P = -ikT (11) 

It is worth mentioning that when Eq.( 11) is satisfied, only the transverse magnetic dipolar scattering will 
be vanishing and consequently the anapole excited (consisting of destructively interfering P and T) can be 
mixed with other excited multipoles, such as magnetic dipoles, electric quadrupoles and so on [16,17,19]. 


3. Results and Discussions 

In Fig. 1 we provide a schematic illustration of anopole excitation in a homogenous particle with plane 
wave incidence. The plane wave is polarized along x direction (in terms of electric field) and propagating 
along z direction. The radius of the particle is R, and the particle is isotropic with refractive index n, or is 
radially anisotropic with transverse (perpendicular to the radial direction) permittivity e* = and radial 
permittivity Sr, with the anisotropy parameter defined as: rj = St/Sr- Both T and P are excited, out of 
phase and of the same magnitude with regard to the scattered power, indicating excitation of the anapole 
[Eqs. (lO)-(ll)]. Also the typical current J and magnetic field H distributions for T only have been shown. 

A. Mixed excitation of anapole and other multipoles within homogeneous isotropic dielectric nanospheres 

At the beginning, we investigate the scattering of homogeneous isotropic {rj = 1) dielectric nanospheres with 
plane wave incidence with R = 200 nm and n = 4 [see the inset of Fig. 2(a)]. The scattered power spectra 
are in Fig. 2(a): the contributions from far field deduced electric and magnetic dipoles P(ai) and M(6i), 
conventional Cartesian electric and toroidal dipoles P and T, and the total scattered power are shown. The 
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Fig. 2: (a) Scattered power spectra for an isotropic (77 = 1) dielectric (n = 4) sphere of iZ = 200 nm (inset). Besides 
the total scattered power (All, black curve), the contributions from current-integrated electric dipole (P, red curve), 
toroidal dipole (T, blue curve), and those from the far-field deduced electric dipole [P(ai), dashed magenta curve] 
and magnetic dipole [M(bi), dashed green curve] are shown. The shaded area is of interest and shown in detail in 
(b). Two points of different A are indicated by black dots (A: A = 1118 nm, B: A = 1075 nm). The corresponding 
near fields on the x — z plane of 1 / = 0 are shown respectively in (A)-(B). The distributions for both the normalized 
Hy (color-plots) and displacement field D = eE on the x — z plane (vector-plots) are shown (the dashed black lines 
denote the boundaries of the particles), as is also the case for Fig. 4. 


shaded region in Fig. 2(a) is the area of interest, and thus is zoomed in and shown in Fig. 2(b), where two 
points A and B are pinpointed and indicated by black dots. At point A of A = 1118 nm, for the transverse 
magnetic dipolar components, there is only toroidal dipole T excitation [P(ai) ss ikT], and the conventional 
Cartesian electric dipole P is completely suppressed; At point B of A = 1075 nm, there is equal out of 
phase T and P excitation in terms of scattering power, leading to an effective anapole excitation defined 
in Eqs. (lO)-(ll) [16,17,19,20]. The corresponding near-held distributions are shown in Fig. 2(A-B): the 
vector-plots correspond to the displacement held D = eE on the x — z plane and the color-plots correspond 
to the normalized magnetic helds along y direction (Hy). We note here that for better visibility we plot 
only transverse magnetic dipolar helds and as is shown: (1) At point A shown in Fig. 2(A), the near held 
exhibits a typical toroidal dipolar distribution, as is schematically shown in Fig. 1 (the magnetic held is 
conhned within the particle and there are accompanying oppositely circulating current loops); (2) Point B 
shown Fig. 2(B) corresponds to an anapole excitation, and at this point the helds outside the particle are 
vanishing, conhrming the negligible overall dipolar scattering and thus effective anapole excitation. 

According to Fig. 2(b), at Point B besides the anapole excitation, there is the coexistence of other multi- 
poles, of which the far held deduced magnetic dipole M(6i) is the most signihcant one. At the same time, 
the scattered power of M(6i) is close to that of both P and T, rendering the anapole excitation of a hybrid 
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Fig. 3: (a) Scattering spectra of radially anisotropic dielectric nanosphere (transverse permittivity is et = 4^ and 
R — 200 nm) showing only contributions from P(ai) and M(bi). The lowest order minimum scattering points of 
P(ai) and M(bi) are pinpointed by circular dots Za and Zb respectively, (b) Spectral dependence of the minimum 
scattering points on rj is shown in (b) and rj = 2.7 is the overlapping point of Za and Zb. 

and mixed nature. Basically to obtain the sole excitation of a pure anapole, the suppression of M(&i) is 
required. 

B. Overlapping the minimum scattering points o/P(ai) and M(hi) within radially anisotropic nanospheres 

It is recently demonstrated that radial anisotropy can be employed to tune spectrally the far field deduced 
electric dipoles P(ai) and thus the anapoles excited [20,28]. In contrast when the transverse permittivity 
Et = is fixed, the spectra of far field deduced magnetic dipoles M(bi) are independent of the anisotropy 
parameter due to its transverse electric nature [28]. In Fig. 3(a) we show the scattering spectra of radially 
anisotropic dielectric nanosphere (transverse permittivity is fixed at Et = 4^ and R = 200 nm). Only 
the contributions from P(ai) and M(bi) are shown. It is clear that P(ai) would blueshift/redshift with 
larger/smaller p, while M(bi) is stable and independent of rj. We indicate the first (counted from large 
wavelength) minimum scattering points of P(ai) and M(bi) by circular dots Za (which corresponds to 
the anapole excitation [20]) and Zb respectively in Fig. 3(a), and their spectra dependence on rj is shown 
Fig. 3(b). According to Fig. 3(b), besides the confirmation of the spectral tunability of P(ai) and stability of 
M(bi), it is further shown that at 77 = 2.7 the minimum scattering points of both P(ai) and M(bi) overlap. 
This means the anapole excitation is achieved with simultaneous magnetic dipole suppression. 
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Fig. 4: (a) Scattered power spectra for an anisotropic (77 = 2.7) sphere of i? = 200 nm and et = 4^ (inset). The total 
scattered power together with the contributions from P, T, P(ai), and M(bi) are shown. Two points C and D of 
A = 907 nm and 852 nm respectively are selected and the corresponding near fields are shown in (C)-(D). 


C. Pure anapole excitation within radially anisotropic nanospheres 

To further show the significance of the overlapping of the minimum scattering points of P(ai) and M(bi) 
and the sole anapole excitation, we show in Fig. 4(a) the spectra at 77 = 2.7 for the radially anisotropic 
nanosphere with et = 4^ and R = 200 nm [see the inset of Fig. 4(a)]. Similar to Fig. 2(a), contributions 
from P(ai) , M( 6 i), P, T together with the total scattered power are shown. Two points C (A = 907 nm) 
and D (A = 852 nm) are indicated by black dots [the corresponding near fields are shown in Fig. 4(C)- 
Fig. 4(D) respectively], whose counterparts (which are similar in terms of both scattering spectra and near- 
held distributions) are respectively points A and B in Fig. 2(a): at point C there is dominant T excitation 
[with magnetic held conhned within the particle and oppositely circulating currents, see Fig. 4(C)] and point 
D corresponds to the excited anapole [with almost zero helds everywhere outside the particle, see Fig. 4(D)]. 
By comparing point B [see Fig. 2] and D [see Fig. 4], the signihcance of introducing radial anisotropy is 
clear: ( 1 ) The anapole is excited at the minimum scattering point of M(bi) and at the same time besides 
the signihcant suppression of M(bi), other multipoles including the electric and magnetic quadrupoles are 
also negligible (not shown here but is clear according to the total scattering), leading to vanishing total 
scattering; (2) Though completely cancel each other, the scattering of P and T is an order of magnitude 
larger than all the other multipoles, showing a very efficient anapole excitation. The two points above are 
sufficient to justify the claim of elusive pure anapole excitation within radially anisotropic nanospheres. 

4. Conclusions and Outlook 

To conclude, in this work we show the elusive excitation of pure anapoles in homogenous nanoparticles, which 
has been realized through employing radial anisotropy of the material permittivity to tune the anapole excited 
to overlap spectrally with the minimum scattering position of magnetic dipoles. At the overlapping position. 
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Cartesian electric and toroidal dipoles can totally cancel the scattering of each other in the far field and at 
the same time the excitation of other multipoles (including but not limited to magnetic dipoles) is negligible, 
proving that the pure anapole has been achieved. 

We emphasize that the anapole excited has been induced by the interferences between Cartesian electric 
and toroidal dipoles, which is featured by null scattering in terms of only transverse magnetic dipolar 
scattering. In a similar way, anapoles characterized by null transverse electric dipolar scattering, or both 
transverse electric and magnetic null scattering of higher orders can be achieved. The investigations into 
scattering nanoparticles from the new perspective of anapole and toroidal multipole excitations can also be 
conducted for other non-spherical structures where numerical methods can be applied to calculate the fields, 
currents, and thus multipoles excited. Also in this work we have confined our study to correcting terms up to 
toroidal multipoles [20,21,24,33]. When the operating spectral region is far beyond the quasi-static regime, 
higher order correcting terms than toroidal multipoles (in terms of kR) will arise inevitably [20,21,24,33], 
and as a result it is expected that anapoles induced by the inferencing between more than two multipoles 
can be obtained. 

Here in this work we show only the demonstration and tuning of anapoles based on radial anisotropy, and 
other effects, such as nonlinear and gain effects, can also be applied for both near and far-held scattering 
manipulations based on anapole control. The investigations conducted into the fundamental Mie scattering 
problem from the new research angle of anapole excitation and tuning is of great signihcance for not only 
deepening our understanding of the basic scattering physics, but also for quite a few scattering related 
applications, such as far-held scattering pattern shaping [28,30,34,35], non-invasive detection and sensing [36] 
and so on. The anapole mode exhibits a number of nontrivial topological properties, which may render new 
opportunities for manipulating electromagnetic waves in an unprecedented way when combined with the 
emerging helds of topological and low dimensional photonics. Moreover, the investigations into anapoles 
should not be conhned to the held of electrodynamics, and similar studies can certainly be conducted in 
other helds such as acoustics, huid dynamics, quantum held theory and so on. 
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